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We report the demonstration of a sensitive absolute grav- 
ity gradiometer based on light-pulse atom interference tech- 
niques. The gradiometer consists of two absolute accelerome- 
ters operated in a differential mode. We report a differential 
acceleration sensitivity of 4 x 10 -9 g/Hz 1//2 and an inferred 
differential acceleration accuracy of less than 10" 9 g. This 
corresponds to a gravity gradient sensitivity of 4 E/Hz 1//2 (f 
E = 10 -9 s -2 ) and an accuracy of better than f E for a 
10 m separation between accelerometers. We demonstrate 
that the instrument can be used to detect nearby masses in a 
vibrationally noisy environment and characterize instrument 
sensitivity to spurious acceleration and rotation noise. 



PACS numbers: 39.20.+q, 03.75.Dg, 04.80.-y, 32.80.Pj 



I. INTRODUCTION 

Precision gravimetry is scientifically and technologi- 
cally relevant. For example, theories which predict vi- 
olations of Einstein's general relativity are manifested 
in gravitational phenomenon such as composition de- 
pendent gravitational forces, time variation of the grav- 
itational constant G, or breakdown of the 1/r 2 law [[!]. 
Some theories predict a fifth force originating from a spin- 
gravity coupling |2j . The gravitational constant itself is 
only known to a few parts in 10~ 4 ||. Technological 
applications lie in the fields of navigation, geodesy, un- 
derground structure detection, and oil and mineral ex- 
ploration. 

A technical problem associated with the characteriza- 
tion of gravitational forces is rooted in the equivalence 
principle: in principle it is not possible to distinguish 
gravitationally-induced accelerations from accelerations 
of the reference frame of measurement. In other words, 
gravimetry is fundamentally challenged by the inability 
to distinguish platform vibrations from gravitational ac- 
celerations. It is well known that this difficulty can be 
overcome with gradiometric measurements in which two 
simultaneous, spatially separated acceleration measure- 
ments are made with respect to a common reference plat- 
form. The difference between these acceleration measure- 
ments suppresses platform vibration noise as a common 
mode, but preserves gravitationally-induced differential 
accelerations (gravity gradients). This non-local mea- 
surement of the curvature of the gravitational potential 
circumvents the equivalence principle and allows for char- 
acterization of gravitational anomalies in vibrationally 
noisy environments. For this reason gravity gradiomc- 
ters are of technological interest. 



In this article we present a sensitive absolute grav- 
ity gradiometer based on atom interference methods. 
We demonstrate a differential acceleration sensitivity of 
4 x 10~ 9 g/Hz 1 / 2 (corresponding to an inferred perfor- 
mance of 2.8 x 10~ 9 g/Hz 1 / 2 per accelerometer). We 
characterize instrument immunity to rotational and vi- 
brational noise through laboratory shake-tests. We mea- 
sure the gravitational gradient induced by nearby mass 
distributions and characterize instrument accuracy. In 
previous work we made a proof-of-principle measurement 
of the gravitational gradient of the Earth using a substan- 
tially less sensitive instrument [Q. 

The central idea behind our atom interferometer grav- 
ity gradiometer is to make acceleration measurements on 
two vertically separated laser cooled ensembles of cesium 
(Cs) atoms in free-fall using a pair of vertically propa- 
gating laser beams. The propagation axes of these laser 
beams are aligned to pass through both ensembles. The 
light-pulse atom interference method is used to measure 
the acceleration of each ensemble with respect to a ref- 
erence frame defined by the phase fronts of the interro- 
gating optical fields. The difference between the mea- 
sured acceleration of each atom ensemble, divided by 
their separation, is a measure of the in-line component 
of the gravity-gradient tensor Ty. (This tensor charac- 
terizes the gravitational field inhomogeneity induced by 
non-uniform mass distributions.) Accelerations of the 
common reference frame - defined by the optical field 
phase fronts - are rejected as a common-mode in the dif- 
ference. 

Each acceleration measurement is accomplished by 
driving an optical pulse tt/2—tt—tt/2 atom interferometer 
sequence between the 6S1/2, F = 3 and F = 4 ground- 
state hyperfine levels of atomic Cs |j| . This method has 
recently been used by Chu and coworkers to measure g 
at the part per billion level, with a sensitivity of 2 x 10~ 8 
g/RV/2 §. 

The interferometer theory is described in detail else- 
where H0]. Here we review the essential results. Each 
interferometer pulse couples two ground-state hyperfine 
levels via a two photon, velocity-sensitive Raman cou- 
pling ||. In the short, intense pulse limit, the opti- 
cal phase associated with the Raman coupling is im- 
printed on the atomic center-of-mass wavepackets. This 
phase locates these wavepackets with respect to the ref- 
erence frame defined by the optical field. The following 
rules, derived from solutions of the Schrodinger equa- 
tion for a two-level atom in the rotating wave approxi- 



mation, govern this imprinting process: |3) 
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ing field at the mean position x of the wavepacket at the 
time t of the interferometer pulse, and states |3) and |4) 
are the two resonantly coupled hyperfine states. Explic- 
itly, <p(t) = k e ff • x(t) + 4>q (t). Here the Raman propaga- 
tion vector k e ff is defined as ki — k.2 ; ki and k2 are the 
wavevectors for two Raman lasers. 

The probability following the three pulse interferom- 
eter sequence for the atoms to be found in |4) (if they 
initially are prepared in |3)) can be evaluated from 
straightforward application of the above rules and is 
given by P(|4)) = [1 - cos A<j)}/2. The phase shift 
A<p = (f)(ti) — 2<p(t 2 ) + 0(^3), where ti is the time of 
the \ th pulse. Evaluating (j){ti) as in [|| gives A</> = 
k e ff • aT 2 + A(/)o . T is the interrogation time between 
interferometer pulses, a is the mean acceleration experi- 
enced by the atoms with respect to the optical fields, and 
A0 O = 4>o (ti) — 2^0(^2) + 4>o (£3)- F° r vibrationally quiet 
environments, 4>o{ti) = ^0(^2) = ^0(^3)- In this case, 
A(f> Q — 0, and measurement of the transition probabil- 
ity following the three pulse sequence determines a. In 
noisier environments, A0 O is no longer zero, and this spu- 
rious shift contaminates an individual acceleration mea- 
surement. 

The gradient signal is derived from two simultaneous, 
spatially-separated, acceleration measurements which are 
made with respect to the same set of Raman laser fields. 
This is accomplished by simultaneously measuring the 
fraction of atoms excited by the pulse sequence at both 
spatial locations. The differential acceleration is given by 
the differential phase shift, Acj) = k e ff • AgT 2 , between 
the upper and lower atom ensembles, where Ag is the 
difference in the mean gravitational acceleration at the 
two accelerometers. In the framework of this analysis, 
low frequency platform accelerations introduce common 
accelerations <5a which cancel in the differencing proce- 
dure. High frequency vibrational noise, which shows up 
in the A</>o terms of each accelerometer, also cancels. 

Our approach differs from that used in other instru- 
ments. First, it uses atoms as proof masses rather than 
macroscopic objects. This eliminates variability from 
device-to-device and provides insensitivity to many envi- 
ronmental perturbations, e.g. temperature fluctuations 
and magnetic fields. The gradiometer references its cal- 
ibration to the wavelength of the measurement laser, 
which is locked to an atomic spectral line, providing abso- 
lute accuracy and long term stability. Because the sensi- 
tive axis is defined by the Raman propagation vector, and 
the acceleration measurements are referenced only to one 
retrorcflector, the two accelerometers may be placed far 
apart without sacrificing common-mode vibration rejec- 
tion. Increasing separation between accelerometers lin- 
early increases the sensitivity to gravity gradients and 
provides insensitivity to near-field perturbations. 

The gradiometer presented here performs favorably 
when compared with existing gravity gradient sensors. 
The state-of-the-art mobile gradiometer is a device that 
uses mechanical mass-spring accelerometers on a rotat- 
ing, gimballed fixture (Lockheed Martin UGM) |J . This 



device has recently been used to perform airborne sur- 
veys of geophysical mass anomalies [[[o) and to charac- 
terize the gravitational anomalies from man-made un- 
derground structures fuJl . Higher sensitivity labora- 
tory devices based on superconducting transducers have 
achieved significantly better sensitivities (< 0.1 E/Hz 1//2 ) 
Jl^pl , However, these devices suffer from tare effects in 
the superconductors (l4[ and their reliance on cryogens. 
More recently, the differential acceleration of two falling 
corner cubes has been measured using a Michelson inter- 
ferometer configuration jn|. However, the demonstrated 
sensitivity of this device of 450 E/Hz 1 / 2 is not competi- 
tive with the Lockheed Martin UGM, and it remains to 
be seen if this device can demonstrate the common-mode 
noise rejection required for mobile applications. 

This remainder of this paper is organized as follows. 
Section II details the apparatus. Section III contains re- 
sults characterizing the instrument's performance. Sec- 
tion IV draws comparisons with related atom interference 
methods. Finally, future work is detailed in Section V. 

II. APPARATUS 

A. Apparatus Overview 

The apparatus is similar to that described in Ref. [Q, 
but several changes have been made in order to achieve 
the current sensitivity. The apparatus consists of two 
laser cooled and trapped sources of Cs atoms. The atoms 
are launched on ballistic trajectories and prepared in a 
particular internal state with optical and microwave tech- 
niques before undergoing the interferometer sequence. 
Following the interferometer sequence, atoms are de- 
tected using a normalizing detection method. Each of 
these stages is described below. In addition, we detail 
the operation of an actively controlled vibration isola- 
tion system which is used to validate our data acquisition 
methodology. 

B. Laser Cooled Atomic Sources 

Each laser cooled atomic source consists of an ultra- 
high vacuum system with a Cs source which maintains a 
low Cs vapor pressure for a magneto-optical trap (MOT) , 
along with the required laser beams and magnetic field 
configuration to form the MOTs. The lasers are delivered 
to the vacuum chambers via optical fibers from a com- 
mon high-power, frequency stable laser system. The laser 
beams are configured in the six beam [1,1,1] geometry (3 
mutually orthogonal pairs) which allows clear access for 
the Raman interferometer beams along the the vertical 
axis. 

In order to obtain good interferometer signal-to-noisc 
ratios (SNRs), it is critical to load atoms quickly into 
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the MOTs in a way which minimizes atom number fluc- 
tuations during the loading process. A grating stabilized 
diode laser (New Focus Vortex) replaces the previously 
used distributed Bragg reflector (DBR) laser as the mas- 
ter laser for the trapping system. The extended cavity 
of the Vortex laser gives it an intrinsic linewidth approx- 
imately ten times smaller than the DBR laser. The Vor- 
tex laser is locked to a Cs transition via standard satu- 
ration spectroscopy techniques in a Cs cell. The locked 
laser has a stability of 1 kHz/Hz 1 / 2 and an instantaneous 
linewidth of 300 kHz. The more stable master laser re- 
duces the shot-to-shot rms number fluctuations of the 
MOTs by a factor of five to a SNR of 200:1. In addi- 
tion, frequency-induced noise during the detection of the 
atoms following the interferometer is substantially sup- 
pressed with the use of this laser system. 

The master laser injection locks two amplifier lasers 
which then seed two 500 mW tapered amplifier lasers 
(SDL 8630E). All of the injection locks transfer the fre- 
quency stability of the master laser, without addition of 
extra frequency or intensity noise. Each tapered amplifier 
laser is split into six equal outputs to provide the twelve 
trapping beams for the two MOTs. The splitting is ac- 
complished using robust and compact free-space fiber op- 
tic beamsplitters (Optics for Research FiberBench). The 
stability of the fiber mounts maintains the splitting ra- 
tio to greater than 1% over many months with no ad- 
justment. The fibers are polarization maintaining with 
an extinction ratio of greater than 30 dB. Polarization 
stability is further maintained with the use of clean-up 
polarizing optics at the output of the fibers (convert- 
ing polarization drift in the fiber to negligible intensity 
fluctuations). The intensity and polarization stability of 
the MOT trapping beams contribute greatly to the re- 
duction of the MOT number fluctuations. The trapping 
beams from the twelve fibers propagate uncollimated to 
the two MOTs. By not collimating the beams, we are 
able to circumvent window apertures limiting the beam 
size. Larger beam waists in the loading region lead to 
higher atom loading rates. The approximate beam waist 
at the trapping position is 2.5 cm (1/e 2 ), and the inten- 
sity is about 1.2 I sot per beam (I sat = 1.1 mW/cm 2 for 
the Cs cooling transition). In this configuration, each 
MOT loads approximately 2 x 10 8 atoms in 1 s. 

After loading the MOTs from a thermal vapor for 1 
s, the cold atoms are launched in ballistic, atomic foun- 
tain trajectories. The launch is accomplished as follows. 
First, the detuning of the MOT beams is switched from 
-10 MHz to -20 MHz from the cooling transition. Next, 
the trapping quadrupole magnetic fields arc turned off, 
and the atoms are held in the -20 MHz detuned op- 
tical molasses for 30 ms while the eddy currents from 
the field switching damp out (the vacuum chambers are 
aluminum). For each MOT, following this holding pe- 
riod, the frequency of the upper three molasses beams is 
ramped down by 1 MHz over 5 ms [by applying the ap- 
propriate RF frequency shifts to an acousto-optic mod- 
ulator (AOM)], while the frequency of the lower three 



beams is ramped up by an equal amount. This frequency 
ramp smoothly transfers the atoms to an optical molasses 
moving vertically at 1.5 m/s. After this ramp, the mean 
frequency of the trapping beams is ramped down to -40 
MHz detuned (still in a moving molasses) over 5 ms, and 
then the intensity is ramped down to half intensity in 
1 ms, held for 0.5 ms, and ramped completely off in an 
additional 0.5 ms. These ramps cool the atoms to 2.3 
fiK (measured with stimulated Raman velocimetry). The 
frequency and intensity ramps are accomplished using 
digitally synthesized waveforms from Hewlett Packard 
HP8770A arbitrary waveform generators (AWGs). 

Following this launching and cooling phase, the cold 
atoms move in a 320 ms, 12 cm high fountain during 
which they are prepared in a special state and then in- 
terrogated by the interferometer sequence. 

C. State Preparation 

Before the atom interference pulse sequence, a state 
selection sequence prepares the atoms in the magnetically 
insensitive my = sublevel and velocity selects an atomic 
ensemble with a velocity spread which is matched to the 
velocity selectivity of the Raman pulse sequence. This 
state selection, which is accomplished with a sequence of 
microwave and optical pulses, is important for obtaining 
high fringe contrast. The details associated with this 
state selection are discussed below. 

Following their launch, atoms are initially distributed 
nearly evenly among the magnetic sublevels of the F 
= 4 ground state. Three orthogonal pairs of magnetic 
field coils, roughly in a Helmholtz configuration, zero the 
Earth's magnetic field and apply a vertical bias field of 
~ 100 mG. This bias allows selective addressing of indi- 
vidual F = 3 to F = 4 ground state hyperfine transitions 
with a 9.2 GHz microwave field (delivered to the atoms 
through a Narda 640 X-band gain horn). First, a mi- 
crowave composite n pulse transfers atoms from the F = 
4, m.f = to the F = 3, my = sublevel. (Composite 
pulse sequences accomplish the population transfer as- 
sociated with ideal ir in an experimentally robust way, 
as described in the subsection below.) Next, a near- 
resonant pulse from the upper trapping beams tuned 
slightly above the 6S!/2, F =4 — > 6P 3 / 2 , F' = 5 transition 
clears atoms in the remaining F = 4 sublevels (via the 
scattering force). Another composite microwave ir pulse 
then returns F = 3, my = atoms to F = 4, my = 0. An 
optical velocity selective composite Raman ir pulse is now 
applied which transfers F = 4, my = atoms within the 
velocity range encompassed by the Raman pulse envelope 
to F = 3, my =0. Finally, a second, near resonant, blue 
detuned pulse clears away the remaining F = 4 atoms. 
At this point, the remaining (state prepared) atoms arc 
ready for use in the interferometer. 

A considerable fraction of atoms are eliminated from 
the initial ensemble in this process: ^8/9 from the in- 
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ternal state selection and another ~ 2/3 in velocity se- 
lection, leaving roughly 4% for the interferometer. 

1. Enhanced Optical Pumping 

In order to increase the m/ = population, an en- 
hanced optical pumping scheme has been implemented. 
This method essentially recycles atoms which are not ini- 
tially in the m/ = target state by a cyclic optical pump- 
ing sequence. In practice we have seen as much as a factor 
of 3 improvement in usable atoms with this method. 

The details of this scheme are as follows. First, a com- 
posite microwave tt pulse is applied to drive atoms from 
the F = 4, to. j = to the F = 3, my = state as be- 
fore. Then a de-pumping beam tuned to the F = 4 to 
F' = 4 transition is applied to optically pump the re- 
maining atoms from F = 4, m.f ^ to F = 3 (where the 
prime still represents the 6P3/2 excited manifold). The 
de-pumping process redistributes the atomic population, 
with approximately 1/7 of the remaining atoms ending 
up in iaif = 0. The process is then reversed with a com- 
posite microwave tt pulse transferring F = 3,m/ = 0to 
F = 4, mj = followed by application of a repumping 
beam to the F = 3 to F' — 4 transition. In principle, this 
entire sequence could be repeated many times, resulting 
in 100% pumping efficiency. 

In practice, inefficiency of the microwave pulses (de- 
scribed below), heating due to spontaneous emission in 
the pumping sequence, and the availability of only a finite 
amount of time to execute the sequence limit the process 
efficiency. We realize a factor of three improvement with 
just one cycle of the sequence described in the previous 
paragraph. 

For vapor cell loaded traps, the overall efficiency is 
also limited by the presence of the background atomic 
vapor. In this case, the de-pumping photons can excite 
atoms in the background vapor which then emit light 
at the repumping frequency. These rescattered photons 
then redistribute atoms in the F = 3, m/ = state to 
other nif 7^ states. Thus there is a trade-off between 
background vapor pressure (which sets the loading rate) 
and the overall efficiency of the scheme (which works best 
at low vapor pressure) . For our operating parameters, we 
typically realize a factor of 1.5 - 2 improvement. 

Finally, we note that this method substantially in- 
creases shot-to-shot atom number fluctuations (from 
200:1 to 60:1). However, our normalized detection 
method is able to effectively suppress the impact of these 
fluctuations on the interferometer signals. 

2. Composite Pulse Techniques 

The above state preparation methods work best with 
efficient coherent population transfer between F = 3 and 
F = 4 states. Less than unit transfer efficiency during 



a standard tt pulse between the ground states can re- 
sult from an inhomogeneous Rabi frequency of the mi- 
crowave or optical pulse seen by the atoms, as well as 
by detunings due to the velocity spread of the atoms. In 
our apparatus, the microwave tt pulses are typically only 
80% efficient due to inhomogeneous field strengths across 
the atom clouds (since we drive the microwave transition 
with horns located outside the vacuum chamber). Fur- 
thermore, the state selection and optical pumping require 
a series of one optical and four microwave tt pulses in the 
two separate chambers. With the current system, it is 
difficult to match the pulse conditions for all pairs of 
pulses due to different microwave intensity gradients in 
each chamber. 

In composite pulse sequences fl6|| , a standard tt pulse is 
replaced with a sequence of pulses with variable area and 
relative phase. In our work, we employ a tt/2 — ttqo° —tt/2 
pulse sequence in place of a tt pulse. The subscript 
90° indicates that the phase of the center tt pulse is 
shifted 90° relative to the tt/2 pulses. The use of this 
sequence increases the transfer efficiency of a pulse for in- 
homogeneous distributions of Rabi frequency and detun- 
ing across the atomic ensemble. Employing these pulses 
for the microwave state preparation pulses increases the 
transfer efficiency from 80% with a regular tt pulse to 
95%. Fig. ^ compares a frequency scan with regular and 
composite tt pulses. More advanced interferometer pulse 
sequences may benefit from the use of composite pulses 
for the optical pulses but are not employed at the mo- 
ment (see Section VI. A). 

In comparison with the adiabatic rapid passage tech- 
nique (ARP) |ri^| , composite pulses are easier to imple- 
ment experimentally, and require less time (or total pulse 
area) to achieve efficient transfer. For example, we find 
that we need approximately 5tt total time (here time is 
referenced to the time to drive a tt pulse at the maxi- 
mum Rabi frequency) for ARP to achieve results similar 
to those achieved with only 2tt time for the composite 
sequence. 

D. Atom Interferometer 

Following its launch and state preparation each atom 
ensemble is subject to the tt/2 — tt — tt/2 interferometer 
pulse sequence. Key experimental details associated with 
this sequence include the frequency stability of the lasers 
used to drive the Raman transition, optimization of the 
Raman excitation parameters, and the physical geometry 
for beam delivery. 

1. Raman Lasers 

As the acceleration-induced phase shifts depend crit- 
ically on the phase and wavevector associated with the 
laser beams used to drive stimulated Raman transitions, 
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it is important to evaluate the possible contributions of 
laser frequency noise on the relative stabilities of the in- 
terferometer signals. We first consider the constraints on 
laser frequency stability, then describe the laser system 
we use to meet these constraints. 

In our excitation geometry, the two ensembles are sep- 
arated by ~ 1.4 m. The Raman fields propagate in an 
asymmetric way to these ensembles (see Fig. |2j). To see 
this, consider the propagation paths of the optical fields 
following the beam dividing optic which separates the 
two Raman fields. The beam of frequency i>i propagates 
roughly x\ 0.3 m before it passes through the first en- 
semble of atoms, while it propagates roughly x\ s» 1.7 
m before it passes through the second ensemble. On the 
other hand, the beam of frequency V2 propagates x\ rs 
3.7 m before it passes through the first ensemble, while 
it propagates x\ ss 2.3 m before it passes through the 
second ensemble. 

If the frequency of the lasers drifts on a time scale short 
compared to the interrogation time T between pulses, 
this can cause an asymmetric phase shift to be read into 
the atomic coherences due to this path asymmetry. For 
example, suppose the laser frequency jitter is Sv, while 
the differential path length travelled by the Raman lasers 
for the two ensembles is £ = (effective path to ensemble 
1) - (effective path to ensemble 2) = (x\ — x\) - {x\ — 
x\) ~ 2.8 m. This leads to a differential phase noise 
of S(f>i aser ~ (k£)(8v/v). For a target interference SNR 
of 1000:1, Sipiaser 1 mrad. For our parameters, this 
implies 8v < 20 kHz. 

In order to achieve frequency stability at this level, a 
second Vortex (external cavity, grating stabilized) laser 
is used as the master laser for the Raman laser system. 
This laser is locked to the 6S 1/2 , F = 3 -> 6P 3/2 , F' 
= 2 crossover resonance (via modulation transfer spec- 
troscopy) using several AOMs to offset the frequency 
to obtain the desired detuning. The lock is maintained 
through the use of a digital signal processor (DSP, Spec- 
trum Signal Indy TMS320C32). The DSP processes the 
lock error signal through a highpass and lowpass channel, 
each operating at a sampling rate of 25 kHz. The high 
and low frequency channels provide feedback to the laser 
current and to the cavity piezo element respectively. Due 
to the presence of long term drift in the piezo element a 
third, very low frequency channel is provided through a 
GPIB command to the laser controller. The measured 
stability of the laser is comparable to that of the mas- 
ter laser used for the optical molasses (and is primarily 
limited by a 5 kHz resonance of the laser's piezo-electric 
transducer). 

The Vortex laser directly injection locks a 150 mW 
(SDL 5422) slave diode laser. This laser is then shifted 
up and down in frequency by 4.756 GHz (160 MHz above 
half the Cs clock frequency) with a high frequency AOM. 
The diffracted orders are then used to injection lock two 
more 150 mW slave laser diodes at a frequency 700 MHz 
red detuned from the F = 3 -> F' = 4 and F = 4 -> F' = 4 
transitions, respectively p8l. The frequency noise on the 



master Raman laser exists on both the slave lasers, but 
their frequency separation remains fixed at the Cs ground 
state hyperfine transition. The detuning from the F' = 
5 level reduces the effect of spontaneous emission due to 
off-resonant single photon excitations from each Raman 
beam. 



2. Raman Beam Delivery 

The Raman beams are delivered to the vacuum cham- 
bers with a polarization maintaining optical fiber in order 
to increase the pointing stability of the Raman beams as 
well as to spatially filter them. Prior to coupling into the 
fiber, the two Raman beams are double-passed through 
80 MHz AOMs. These AOMs are controlled by another 
HP8770 AWG which allows dynamic frequency, phase, 
and intensity tuning of the Raman beams during the in- 
terferometer pulse sequence. During the interferometer, 
the frequency of the Raman lasers must be chirped in 
order to maintain a resonance condition with the accel- 
erating (Doppler shifting) atoms. 

The two Raman beams are overlapped with orthogonal 
linear polarizations on a polarizing beamsplitting cube 
and passed through a Pockel's cell polarization modula- 
tor (ConOptics 350-50) after the double-passed AOMs. 
The Pockel's cell is used to reverse the effective Raman 
laser propagation direction, as described below. The Ra- 
man beams then are coupled into a polarization main- 
taining fiber with 75% efficiency and sent to the gra- 
diometer. For wavefront quality, after the fiber the Ra- 
man beams are collimated with a 1.1 cm focal length 
aspheric lens in conjunction with a high surface quality 
50 cm focal length spherical lens. This lens combination 
results in a uniform phase front for the Raman beams. 
After collimation the Raman beams have a 1.0 cm (1/e) 
beam waist. All optics in the Raman beams' propaga- 
tion path after the optical fiber are of high surface figure 
(A/10 or better) in order to preserve the phase front ho- 
mogeneity of the Raman beams. 

After collimation, the Raman beams enter a racetrack 
geometry in order to obtain counterpropagating beams 
for the Raman interaction (see Fig. g). The racetrack 
configuration starts with a polarizing beamsplitter cube 
that separates the two orthogonally polarized Raman 
beams. The two Raman beams then parallel propagate 
vertically through the vacuum chambers with one beam 
passing through the axis of the atom ensembles and the 
other 2 cm off-axis. After passing through both cham- 
bers, a corner cube retroreflector redirects the off-axis 
Raman beam to counterpropagate with the on-axis Ra- 
man beam, resulting in two counterpropagating beams. 
The use of the corner cube decreases the tilt sensitivity 
of the apparatus by keeping the Raman beam propaga- 
tion axis constant as the cube is subjected to spurious 
tilts. In this racetrack, standing waves are eliminated, 
which maintains overall intensity stability by suppressing 
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etalon effects from the Raman beams. Also, spontaneous 
emission is reduced by half compared to using collinear 
retroreflected beams. 



3. Raman Propagation Reversal 

A technique to reduce many systematic interferometer 
phase shifts involves reversing the effective Raman prop- 
agation vector k e ff = ki — k>2 . Because the gravitational 
phase shift is proportional to k e ff • g, reversing the sign 
of k e ff changes the sign of the gravitational phase shift. 
However, several systematic phase shifts, such as second- 
order Zeeman shifts from magnetic fields and any resid- 
ual AC Stark shifts, have no dependence on the Raman 
wavevector direction. Subtracting the phases obtained 
from consecutive experimental cycles using two reversed 
propagation directions gives twice the gravitational phase 
shift, but removes these systematic shifts. The propa- 
gation reversal is accomplished using the Pockel's cell. 
The Pockel's cell rotates the polarization of both Raman 
beams by 90° when activated. This rotation causes the 
direction the Raman beams take through the racetrack 
to switch, i.e. k e ff — > — k e ff . Propagation direction can 
be switched from one shot to the next. 



4- Raman Beam Parameters 

A theoretical study of the interferometer contrast as 
a function of Raman beam waist and detuning was per- 
formed. A finite Raman beam size gives rise to a spa- 
tially inhomogeneous Rabi frequency across the atom 
cloud, causing dephasing during the interferometer. Sim- 
ilarly, the velocity spread of the atoms along the Raman 
beams causes inhomogeneous broadening due to differ- 
ential Doppler shifts across the atom ensemble. Typical 
Rabi frequencies are around 30 kHz, and the initial width 
of the thermal spread of the atom ensemble is about 45 
kHz. Incorporating the finite beam size and the initial 
atomic velocity distribution function into our analysis, we 
find an optimal 35% contrast for a Raman beam waist of 
1.0 cm radius (1/e), and a detuning of » 1 GHz from F' 
= 5 (our chosen operating point). 

In addition to the interferometer contrast, systematic 
AC Stark shifts were studied. AC Stark shifts from the 
Raman pulses themselves cause spurious phase shifts in 
the interferometer if unconstrained. However, with a 
two-photon Raman transition, the AC Stark shift is the 
difference between the individual AC Stark shifts from 
each beam and can be zeroed by adjusting the ratio be- 
tween the two Raman beams. The Stark shift is balanced 
with a beam intensity ratio of ~ 1.6:1 for the chosen 
Raman detuning. The Stark shift is balanced empir- 
ically by inserting off-resonant Raman pulses within a 
microwave tt/2 — n — tt/2 interferometer and adjusting 
the Raman beam intensity balance to zero the optically 



induced phase shift. This ratio of 1.6:1 for the chosen Ra- 
man detuning agrees well with theoretical predictions. 

5. Interferometer Operation 

The gradiometer is typically operated in its most sen- 
sitive configuration with the interferometer pulses at a 
maximal spacing of T = 157.5 ms. This time is limited by 
the vacuum chamber size which constrains the fountain 
height to 12 cm. Following the three-pulse interferom- 
eter sequence, the population distribution of the atoms 
in each ensemble is measured. In order to extract the 
gravitationally- induced phase shift, the phase of the final 
interferometer pulse needs to be scanned. The scanning 
is accomplished digitally with the AWG by scanning the 
phase of the RF waveform applied to the low frequency 
AOMs in the Raman beam paths. This phase scanning 
is in addition to the frequency chirp applied during the 
pulse. The total cycle time is 0.7 - 1.4 s, depending on 
the trap loading parameters. 

E. Detection System 

1. Detection Apparatus 

In order to realize high interferometer sensitivity, 
atoms must be detected with a high SNR following the in- 
terferometer pulse sequence. Our detection system uses a 
balanced, modulation transfer technique to reduce laser- 
induced detection noise and differentiate cold atoms from 
thermal background atoms. This detection method is de- 
scribed in [|l9| and is summarized here. 

Balanced detection uses two parallel, horizontally 
propagating probe beams, separated vertically by 2.5 cm, 
to detect atoms in one or both states following the inter- 
ferometer. The absorption of the probe beams by the 
atoms is detected on a balanced photodiode with two in- 
puts that are subtracted before amplification. Typically, 
the absorption by the cold atoms is around 0.1%. 

We use two detection pulses as follows. Following the 
interferometer, atoms in the F = 4 state are detected by 
pulsing on near-resonant probe beams for 4.6 ms when 
the atom cloud is in the upper detection beam. During 
this pulse, the mean velocity of the F = 4 atoms is slowed 
to rest by optical forces induced by the detection beams. 
Atoms in the F = 3 state continue to fall, however. When 
these atoms intersect the lower probe region, atoms in the 
F = 4 state have not moved out of the upper detection re- 
gion. Pulsing on the probe beams for another 4.6 ms fol- 
lowing a short repumper pulse measures the differential 
absorption between the two hyperfine states. These two 
pulses provide sufficient information to infer the overall 
transition probability in a way which is immune to shot- 
to-shot fluctuations in the total atom number. We have 
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shown that this method is immune to laser frequency and 
amplitude technical noise as well Jll|. 

In addition, a modulation transfer technique is em- 
ployed to remove noise from the thermal background va- 
por. Here a frequency modulated (FM) pump beam is 
applied along a vertical axis orthogonally to the probe 
beams during each detection pulse. The nonlinear in- 
teraction between the FM pump and the atoms modu- 
lates the complex index of refraction of the atoms. This 
modulation produces an amplitude modulation (AM) on 
the probe beams pfj| , pT| . The AM is detected by the 
balanced detector, mixed down to DC with a double- 
balanced mixer, and integrated with a Hewlett Packard 
3458A multimeter. The orthogonality of the pump and 
probe beams provides a velocity selectivity that allows 
rejection of signals from the fast-moving thermal back- 
ground atoms. 

2. Detection Noise Analysis 

The amplitude of the modulated absorption signal at 
the photodiode is about 0.8 pW per atom. The detection 
photodiodes have noise-equivalent powers corresponding 
to a 60 atom detection sensitivity in the 4.6 ms measure- 
ment window. The digital voltmeters that acquire the 
signal from the photodiodes are only slightly noisier with 
a noise floor corresponding to 100 atoms. The largest 
contribution to the intrinsic noise comes from the detec- 
tion light itself. Because absorption detection is used and 
the absorption is small, there is a substantial amount of 
unabsorbed light, > 99%, striking the detector. Shot 
noise on the number of photons incident on the detec- 
tor during the integration window is the leading intrinsic 
noise source. The shot noise power is 0.25 nW, resulting 
in a minimum detectable signal of ~ 300 atoms. This 
noise dominates the technical noise sources. 

There is also a noise component of similar size due 
to a small number of atoms in undesired states which 
survive the state preparation. Slight changes in laser 
frequency and selection pulse efficiency during the state 
preparation cause this number to fluctuate. However, 
this noise source is common between the two chambers 
and is also suppressed by the balanced detection method. 

To summarize, the noise of our detection system allows 
detection of transition probability at the atom shot-noise 
limit when we detect more than ~ 10 5 atoms. These 
noise limits are discussed in the context of our overall 
instrument sensitivity below. 

F. Signal Extraction 

1. Interference Fringe Fitting 

As described above, the phase of the final interferom- 
eter pulse is scanned electronically. A straightforward 



method of extracting gravity gradient information is to 
determine the gravitationally induced phase shifts in each 
atom ensemble by performing least squares sinusoidal 
fits on the observed interference fringes. This is possi- 
ble when vibration-induced phase noise is < 1 rad. The 
gravity gradient is obtained by subtracting the two phase 
shifts from each other. Vibrational phase noise and lo- 
cal oscillator phase noise cause the phase extracted by 
the sine fits to be shifted. However, these noise sources 
couple to the two accelerometers in a common-mode way. 
This common-mode behavior results in the two sinusoidal 
fits being shifted by an identical amount, and any effect 
of common-mode noise is cancelled in the subtraction 
used to obtain the gravity gradient. We study the statis- 
tics of the resulting phase differences under static gravity 
gradient conditions to estimate instrument noise. 

We find that the distribution of the residual noise con- 
tains outlying points. Eliminating these outlying points 
increases the SNR by up to a factor of six. The reduc- 
tion of the number of points is incorporated into the data 
collection time in determining the sensitivity. We are 
presently studying possible sources of this noise compo- 
nent. 

The ratio of the interference fringe amplitude to the the 
standard deviation of the phase difference distributions 
determines the instrument SNR. The side of a fringe, 
i.e. the linear slope of a sine wave, is most sensitive to 
phase shifts, with a sensitivity given by 8(f) = 1/SNR. 
For gravitationally induced phase shifts, the sensitiv- 
ity to a change in the gravitational acceleration is Sg 
— <50/(2kT 2 )t 1 / 2 , where t is the data acquisition time 
to achieve the uncertainty S(f>. Dividing by the chamber 
separation determines the sensitivity to gradients. 

2. Magnetic Phase Shifting 

Good common-mode noise suppression requires that 
the lower and upper chamber fringes be acquired in 
phase. However, the Earth's gravity gradient of ~ 3000 
E will cause a relative phase shift of ~ 1.5 rad between 
the two chambers. In order to accommodate this shift, 
a bias magnetic field is pulsed on for 67 ms in the lower 
chamber during the atom interferometer. This field pulse 
causes a phase shift due to the second-order Zeeman ef- 
fect. The amplitude of this pulse is chosen to produce 
a shift which compensates the shift due to the gravity 
gradient, allowing both fringes to be acquired in phase. 

We have quantitatively studied the effectiveness of 
common-mode noise suppression as a function of the 
phase alignment of the interferometer fringes. Fig. ^ 
shows a comparison of the predicted dependence of the 
noise on the relative phase between fringes with experi- 
mental SNR data. For simplicity, we scanned the relative 
phase between pulses simply by changing the interroga- 
tion time T. (For this study, we disabled the magnetic 
phase shifter and worked with relatively short interroga- 
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tion times.) The data shown in Fig. || is the SNR of the 
phase difference between the two fringes obtained by per- 
forming least squares fits as described above on a set of 
interference fringes. The prediction is the simulated SNR 
of the phase difference between two phase-mismatched 
sine functions with phase noise added. The predicted 
SNR is added in quadrature with uncorrelated detection 
noise noise in order to compare the simulation with ex- 
perimental data in this curve. As shown, the data is in 
good agreement with the theory, demonstrating the need 
for the magnetic phase shifting pulse. 

The stability of the bias field used in the magnetic 
shifter has also been studied. Noise on this bias field 
could lead to extra phase noise on the gravity gradient 
signals. This noise has been investigated by applying 
a bias pulse in a microwave clock 7r/2-7r/2 experiment 
and comparing the SNR for a weak pulse (phase shift ~ 
1.2 rad) and a strong field pulse (phase shift ~ 67 rad). 
For the large phase shift, no extra phase noise or drift is 
seen at a fractional uncertainty in the phase of less than 

3 x 10 ~ 6 . This implies that for a 1 rad shift, the bias 
pulse contributes noise at a level below 8 x 10~ 13 g. 

3. High Phase Noise Regimes 

In the case where phase noise is greater than 1 rad, 
the noise renders it impossible to characterize instrument 
noise using the least squares method described above. In 
the most sensitive modes of operation, the vibrations of 
the reference platform induce phase noise much larger 
than this level. A different analysis technique must be 
performed in this regime, using a point-by-point anal- 
ysis. After collection of the gradiometer data, a signal 
extraction algorithm removes amplitude and phase noise 
from the signals: Our detection method allows measure- 
ment of the number of atoms in the F = 4 state and 
the population difference between the two states. With 
this information, common amplitude fluctuations in each 
chamber, primarily from number fluctuations from the 
MOT loading, can be removed. After removing ampli- 
tude noise, a cross-chamber normalization is performed 
to reject phase noise that is primarily vibration-induced. 

The details of this noise analysis procedure are as fol- 
lows. After each interferometer cycle (which represents 
one gradient measurement) two samples are acquired in 
each accelerometer: signals proportional to the F = 4 
population and to the differential population (i.e. pro- 
portional to the number of atoms in F = 4 minus the 
number of atoms in F = 3). The two samples are com- 
bined to infer the total number of atoms present in the 
interferometer during each experimental cycle. The F = 

4 signal is then divided by this total atom number to 
remove any fluctuations in the amplitude of the F = 4 
signal from shot-to-shot atom number fluctuations. We 
normalize each interferometer with this procedure. To re- 
move common phase noise between the two chambers, a 



series of experimental cycles is taken, and a least squares 
minimization (via Gaussian elimination) is performed on 
the quantity (Si— aS2 — /3) 2 where Si and S2 are the shot- 
by-shot normalized F = 4 population levels from the two 
interferometers. The fit constants a and (3 are used to 
compensate for possible differences in interference con- 
trast between the two accelerometers. 

The residuals of the Gaussian elimination procedure 
are used to estimate instrument noise, as show in Fig. 
^. Again, this distribution is non-Gaussian, and outly- 
ing points are discarded to obtain SNR and short-term 
sensitivity estimates. 

Note that while this method is suitable for characteri- 
zation of the instrument noise floor under static gradient 
conditions, further work is need to demonstrate effec- 
tive algorithms for extraction of dynamic gradient sig- 
nals (such as would be present in moving-platform appli- 
cations). We are currently exploring algorithms for this 
purpose. 

G. Vibration Isolation Subsystem 

At the most sensitive gradiometer operation, vibration 
phase noise is large, and the high phase noise algorithm is 
employed. However, in order to verify the validity of this 
algorithm, a vibration isolation system was constructed 
to remove most of the vibration- induced phase noise from 
the interference fringes. With this reduction in phase 
noise, the least squares fit algorithms also can be used to 
reduce the data. 



1. Mechanical Design 

The primary object in the instrument that must be 
isolated from vibrations is the Raman beam corner cube. 
All other optics are positioned so that any vibrations 
Doppler shift the two Raman beams in a common way, 
and the Raman difference frequency remains unchanged. 
The Raman beam corner cube retroreflector is mounted 
on a Newport sub-Hertz platform (SHP) which provides 
the principal vibration isolation (see Fig. ||). The SHP 
is guided by a linear air bearing (New Way S4010002) 
along the vertical axis. The SHP provides isolation in 
the range of 0.5 Hz - 40 Hz. An accelerometer (Tele- 
dyne Geotech S-510) is mounted on the SHP to monitor 
platform accelerations. The corner cube is attached to 
the platform by a stack of two pieces of 1 in. thick lead 
filled acoustic foam separated by a 0.5 in. thick sheet of 
aluminum. The double stack of acoustic foam reduces 
vibrations of 30 Hz and higher by more than 20 dB. A 
linear voice coil actuator provides active feedback to the 
SHP. Additionally, the voice coil can be used to drive 
platform accelerations; this shake-testing is the subject 
of Section III.B. 
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2. DSP Servo System 

Here we describe the active servo system for the SHP 
platform. The active feedback loop starts with the ac- 
celcromctcr to monitor vibrations. The accelerometer 
output is processed by a DSP (Spectrum Signal Indy 
TMS320C32), which we use to filter digitally the ac- 
celerometer input (as described below) and generate the 
feedback error signal. The feedback signal, after being 
buffered by a voltage amplifier, closes the feedback loop 
by driving the voice coil mounted between the SHP and 
the platform support. We apply the following digital fil- 
ters in processing the accelerometer signal. First, a lag 
filter with a bandwidth of 1 Hz to 80 Hz rolls off the 
feedback below the accelerometer 's 100 Hz high frequency 
cut-off. Next, a second lag filter with identical bandwidth 
is used to make the gain roll-off second order. Finally, two 
lead filters are applied to keep the system from oscillat- 
ing at low frequency near the closed-loop SHP resonance 
of 0.03 Hz, which is also close to the internal highpass 
frequency of the accelerometer. The two lead filters have 
bandwidths of 38 mHz - 200 Hz and 380 mHz - 200 Hz re- 
spectively. The total gain of all four filters is 1600. This 
work is similar in concept to that reported in Ref. |0 . 

Using this servo, we are able to reduce the vibrations to 
near the noise floor of the accelerometer (10 -8 g/Hz 1 / 2 ) 
over a bandwidth of 40 mHz - 25 Hz. Higher frequencies 
are passively attenuated by the acoustic foam. With the 
addition of the vibration isolation system, phase noise 
from accelerations of the corner cube is reduced to less 
the 1 rad, and least squares sinusoidal fits may be per- 
formed on the fringes for the longest interrogation times. 

H. Microwave Generation 

The generation and delivery of the 9.2 GHz microwave 
field is briefly described here. The microwave field is cou- 
pled to the atoms through RF horns attached to view- 
ports on the MOT chambers. The microwave frequency 
is tied to a 10 MHz reference, temperature stabilized, 
master crystal oscillator (Oscilloquartz OCXO, stability 
of 1.4x10" 13 in 1 s). The reference oscillator drives a 
100 MHz phase locked oscillator (PLO, Wenzel 500-0732) 
which is the input to a Microlambda (MLPE 1162) 9.2 
GHz PLO. The RF is mixed in a single sideband mixer 
with an ~ 7.4 MHz signal from an AWG which is also 
phase locked to the reference oscillator. The AWG is 
used to scan the RF frequency and phase. The mixer 
output is amplified up to ~ 1 W and sent to the horns. 
The RF power is controlled by the AWG output, and the 
relative power to the two chambers can be adjusted with 
appropriate attenuation in the two paths. 

We perform a microwave 7r/2 - 7r/2 clock experiment as 
a diagnostic to check the phase noise performance of our 
oscillators and the noise performance of our detection sys- 
tem. We have shown that we can detect microwave clock 



fringes with 1000:1 SNR using our normalized detection 
|l9| . This SNR is at the atom shot noise limit for our 
fountain with no velocity selection (~ 10 6 atoms/shot). 

III. INSTRUMENT PERFORMANCE 
A. Sensitivity 

Each interference fringe is typically recorded with 15 
consecutive cycles of the experiment. The number of 
points per scan is kept small to decrease sensitivity to 
long term drifts in signal amplitude and contrast. The 
source of such drifts include fluctuations in laser power 
and drifts in the Cs vapor pressure in the vacuum cham- 
bers. The normalized data reduction method described 
in Section II. F is performed on the data. The result- 
ing observed SNR is typically 150:1. This corresponds 
to a differential accelerometer performance of 4 x 10~ 8 
m/s 2 /Hz 1 / 2 , or ~ 4 x 10~ 9 g/Hz 1 / 2 . Scaling this to a 
10 m baseline gives an inferred gradient sensitivity of 4 
E/Hz 1 / 2 . Fig. |I] shows typical high sensitivity fringes in 
the upper and lower chamber for a T = 157.5 ms inter- 
ferometer. The fringe contrast is typically w 33%. 

Fig. |^ compares interference phase scans with the ref- 
erence platform rigidly attached to the optical table 
(~ 10~ 5 g technical acceleration noise) with the refer- 
ence platform servoed using the vibration isolation sys- 
tem. We can use the data acquired under servoed con- 
ditions {e.g. Fig. |TJb) to provide an independent check 
on the observed SNR, by performing the least-squares 
fitting analysis detailed in section II.F.l, and comparing 
it with the point-by-point analysis of II. F. 3 for data ac- 
quired without the servo {e.g. Fig. |^a). We find that 
both approaches yield consistent noise estimates. 

1. Noise 

The dominant noise source is atom shot noise: the 
Poissonian fluctuations that arise from detecting atoms 
in coherent su per position states, also called quantum pro- 
jection noise f23[| . Atom shot noise scales as the square 
root of the number of atoms detected. As long as there is 
a sufficiently large number of atoms contributing to the 
interference fringes, this noise source will be the domi- 
nant noise source. 

For the T = 157.5 ms gradiometer, the demonstrated 
SNR is predominantly limited by atom number shot 
noise. The fringe amplitude corresponds to about 2.5 x 
10 5 atoms. At 33% fringe contrast, there is a mean offset 
of 7 x 10 5 atoms. Based on atom shot noise, these atom 
numbers put a limit of 300:1 on the SNR for each interfer- 
ence fringe (2.5x 10 5 /(7x 10 5 ) 1 / 2 ). Additionally, normal- 
izing each chamber reduces the SNR for each gravimeter 
to (300/2 1 / 2 ):!. Subtracting the two gravimeter signals 







to produce a gradient signal decreases the SNR by an- 
other factor of 2 1 / 2 , fundamentally limiting the SNR for 
the full gravity gradiometer to 150:1. 

The noise produced from photon shot noise on the 
unabsorbed portion of the probe detection light, corre- 
sponding to ^300 atom minimum detectable signal, is ap- 
proximately at the 800:1 level for each individual fringe, 
giving about a 400:1 limit after the normalization and 
subtraction of the two fringes. This means that photon 
shot noise does not impose a substantial SNR limit for 
the current number of atoms, but it would become sig- 
nificant for reduced atom numbers. 



2. Mass Detection 

In order to demonstrate further the sensitivity of the 
gravity gradiometer, measurement of the gradient of a 
nearby object was performed. Previously, the Earth's 
gravitational gradient was measured With the im- 
proved sensitivity, a measurement of the gravity gradi- 
ent from small test masses has been made. For this 
measurement, eight lead bricks (~ 12.5 kg each) were 
stacked symmetrically about the lower chamber, ~ 0.2 
m from the apex of the atomic fountain. The calcu- 
lated acceleration from this configuration of mass should 
yield a8.2xl0~ 9 g signal. The measured acceleration is 
8.1 ± 2.1 x 10~ 9 g which agrees well with the expected 
phase shift. 

We are currently extending this work to measure the 
gradient of a well characterized mass. Our goal is to 
measure G, the gravitational constant, at the part per 
thousand level. 



B. Immunity to Environmental Noise 

An important feature of the gradiometer is the ability 
to reject common mode accelerations in the two measure- 
ments. This is critical for a device which might be used 
on a moving platform. In order to demonstrate this capa- 
bility, we performed a series of acceleration and tilt tests 
to characterize the instrument's sensitivity to platform 
vibration and tilt noise. We characterize the effects of 
accelerations and tilts by measuring instrument SNR as 
a function of the amplitude and frequency of an external 
platform drive. 

1. Linear Acceleration 

We characterized immunity to linear accelerations by 
shaking the platform on which the retroreflecting corner 
cube is mounted. As in earlier work, this corner cube 
plays the role of the acceleration reference for the ac- 
celeration measurements We drive the platform 
by applying sinusoidal drive currents to the voice coil 



coupled to the vibration isolation platform. We monitor 
the resulting platform acceleration with the platform ac- 
celerometer (described in Section II.G.l). For this work, 
the high frequency lead-foam passive isolation between 
the platform and the corner cube was removed. This 
study was done with the interferometer operating in its 
most sensitive configuration (at T = 157.5 ms interroga- 
tion time). Fig. ^ summarizes the results of these tests. 
For frequencies in the range 1 to 100 Hz no significant 
degradation of the SNR was observed for drive ampli- 
tudes up to 2.5 x 10 -2 g. This corresponds to a common- 
mode rejection ratio of 140 dB. At amplitudes greater 
than ~ 0.1 g, the accelerations are large enough that 
individual Raman pulses are Doppler-shifted out of the 
Raman resonance condition. At this point, we observe 
poor interference fringe contrast. 

2. Rotation 

Tilts of the Raman interferometer beams are expected 
to degrade the sensitivity of the measurement. In or- 
der to study the effects of small tilt displacements, we 
floated the optical table on which the gradiometer ap- 
paratus was mounted using commercial pneumatic legs 
(Newport 1-2000). We drove tilt motions with an appro- 
priately placed voice coil, and characterized these mo- 
tions using a tilt meter (Applied Geomechanics 755-1129; 
specified resolution of 1 /irad, 20 Hz bandwidth). The 
rotation vector associated with the tilt motion was in 
a plane perpendicular to the Raman propagation axis. 
Fig. H shows observed SNR versus maximum tilt ampli- 
tude for different driving frequencies. 

Both centrifugal and Coriolis forces lead to a degrada- 
tion in instrument sensitivity. A rotation f2 generates a 
centrifugal acceleration R fi 2 , where R is the distance of 
one accelerometer from the center of rotation. For two 
accelerometers separated by distance SR, the differential 
acceleration is SR fl 2 , while the acceleration gradient is 
fl 2 . For our operating parameters, we expect this mech- 
anism to begin to degrade the SNR at ~ 10~ 4 rad/s. 
However, full loss in SNR is not expected until £1 ~ 10~ 3 
rad/s for this mechanism, much higher than the rotation 
rates employed in this study. The Coriolis force leads to a 
loss in fringe contrast. In this case, the source transverse 
velocity spread of <5v ~3 cm/s leads to an interference 
phase spread of S<j) ~ 2k e// 8v ft T 2 . At O ~ 10" 4 , 5<j) ~ 
1 rad, and we expect a full loss in contrast. This is in 
reasonable agreement with our observations. 

C. Accuracy Estimation 

To demonstrate accelerometer (hence gradiometer) ac- 
curacy, we monitored accelerometer outputs over ex- 
tended periods of time. As in || we observed the daily 
fluctuations in the measured acceleration due to gravity 
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induced by the motion of the Sun and Moon. Fig. || 
shows data from one accelerometer taken over a period 
of several days compared to the Tamura 87 tide model 
pJJ. The only free parameter in this fit is an overall 
phase offset. From this data set, we constrain the ac- 
celerometer phase offset to better than 10~ 9 g over the 
two day measurement cycle. Subsequent measurements, 
taken several weeks later, resulted in a similar determi- 
nation of the fitted offset. From these measurements we 
infer an accelerometer accuracy of < 10 -9 g over time pe- 
riods of days. For a 10 m baseline gradient instrument, 
this corresponds to an inferred accuracy of better than 1 
E. Further accuracy studies are underway. 

IV. DISCUSSION 

In this section we discuss the light-pulse method in 
the context of other de Broglie wave interference meth- 
ods based on laser-atom interactions, as well as alternate 
geometries for gradient measurements. 

A. Interferometer Comparisons 

Here we compare de Broglie- wave gravimeters based on 
multiple pulse techniques using stimulated Raman tran- 
sitions, Bragg scattering, diffraction in the Raman-Nath 
regime, adiabatic transfer, and the AC Josephson effect 
with Bose-Einstein condensates (BEC). 

We have recently demonstrated large-area, multiple 
Raman pulse techniques that may enable higher sensi- 
tivity gravimeters and gradiometers ]25| ] . By using extra 
stimulated Raman transition pulses of alternating propa- 
gation direction, a large relative momentum can be trans- 
ferred to the two interfering atomic wavepackets, result- 
ing in a large-area interferometer. The sensitivity of this 
device to gravitational phase shifts increases linearly with 
the relative momentum imparted to the wavepackets, or 
the number of extra pulse applied. For instance, a gra- 
diometer with one extra set of pulses has a relative split- 
ting of 6fik and a sensitivity to gravity gradients given 
by A</> = 6k • AgT 2 . In principle, each subsequent Ra- 
man pulse should not affect the fringe contrast, except 
for a small amount of spontaneous emission. This process 
should be extendable to a large number of extra pulses 
and a large relative momentum. 

This method might be extended using special se- 
quences of composite pulses (see Section II.C.2). The 
benefits of using composite pulses in this context would 
be an extremely large input velocity acceptance for a 
fixed Rabi frequency and robust suppression of spatial 
inhomogeneities in Rabi frequency (which arise due to 
the finite size of the Raman laser beams). For exam- 
ple, consider the following sequence. A high efficiency 
composite microwave pulse could be used for the initial 
ir/2 pulse. Next a composite optical tt pulse gives the 



wavepackets 4Kk momentum splitting, and after a long 
drift time two composite 7r pulses redirect the wavepack- 
ets back towards each other. Finally, a last composite 
optical pulse and another composite microwave pulse re- 
combine the wavepackets to complete the interferometer. 
Simulations show that this larger area interferometer has 
the potential to have a very high fringe contrast (because 
of the symmetric composite pulses). 

Bragg scattering-based interferometers diffract atoms 
from standing waves of laser light. As in the light-pulse 
interferometers, these interactions can be configured as 
atom optic beamsplitters and mirrors [p6|| . The primary 
virtue of Bragg scattering interferometers is that the 
atoms always remain in the same internal state. This 
reduces the sensitivity to systematic phase shifts such as 
Zeeman shifts and AC Stark shifts, since the wavepack- 
ets in each arm of the interferometer experience symmet- 
ric phase shifts due to these effects. In comparison, the 
Raman method requires the use of propagation vector 
reversal techniques to gain immunity to these possible 
systematics. Furthermore, high-order Bragg processes 
can be used to create large-area, high-sensitivity instru- 
ments. However, high-order Bragg processes operate ef- 
ficiently over a relatively narrow range of initial atomic 
velocities (significantly less than a photon recoil velocity) . 
This severely constrains the fraction of atoms which can 
contribute to the interference signal, and thus the atom 
counting rate. On the other hand, BEC or atom laser 
sources may eventually produce extremely bright atomic 
beams (having excellent velocity collimation) . In this 
case, Bragg processes may become a competitive choice 
for interferometer sensors. 

A gravimeter based on diffraction in the Raman-Nath 
regime has been demonstrated [p7| . In this interferome- 
ter, short, intense pulses of light are applied to the atomic 
ensemble. Due to the shortness and intensity of these 
pulses, atoms are diffracted into a large number of diffrac- 
tion orders. Two such Raman-Nath pulses are used to 
construct the interferometer, and because of the wide 
spread of wavepacket momenta, many different interfer- 
ing paths exist. Some of these paths overlap and interfere 
after an echo time, and the echo time determines the sen- 
sitivity to gravitational phase shifts. The Raman-Nath 
diffraction also populates many higher lying momentum 
states that do not contribute to the closed interferom- 
eter paths, which decreases the interferometer contrast 
and places severe constraints on the SNR. High sensitiv- 
ity gravimeters, which require good SNR for their oper- 
ation, have yet to be demonstrated using this approach. 

Adiabatic transfer recently has been used in a proof- 
of-principlc demonstration of a possible large-area inter- 
ferometer |p9| . In this approach, atoms are put into co- 
herent superpositions of two states using a microwave 
pulse. Momentum is transferred to one state in this su- 
perposition by adiabatically transferring it from the m^- 
= sublevel to the highest (or lowest) lying m/ level. 
Adiabatic transfer is used subsequently to manipulate 
and ultimately to overlap the wavepackets. A final mi- 



ll 



crowave pulse is then used to interfere these wavepackets. 
This method can transfer up to 2mfKk momentum to one 
arm of the interferometer. The primary limitation to the 
utility of this method is that the adiabatic transfer pop- 
ulates magnetic field sensitive sublevels. This makes the 
interferometer sensitive to Zeeman shifts. While such an 
interferometer may demonstrate a high sensitivity, it re- 
mains to be seen whether it can achieve high accuracy. 

Finally, the AC Josephson effect in arrays of Bose- 
Einstein condensed atoms |^8| has recently been used 
to make a proof-of- principle gravity measurement. In 
this approach, condensate atoms tunnel from an array 
of vertically spaced lattice sites. Atoms tunnelling from 
different sites subsequently interfere. The resulting in- 
terference pattern is a periodic train of atom pulses 
whose frequency depends on the strength of the gravi- 
tational potential. This frequency can be measured with 
high accuracy using conventional atom detection tech- 
niques. A major technological drawback to this technique 
is the need for Bose-Einstein condensed atomic sources, 
which still are difficult to produce and not suitable for a 
portable apparatus. Also, the time required to condense 
atoms is typically in excess of 30 s, which results in a low 
instrument bandwidth. However, if techniques to pro- 
duce robust, BEC atom sources improve, this method 
may become viable for future instruments. 

B. Direct Gradient Measurements 

1. Multi-loop Interferometers 

There are several Raman-pulse based schemes which 
can be used for direct gravity gradient measurements. 
The simplest of these is the double-loop, or figure-eight, 
interferometer |p0[| . This geometry can be achieved in 
the light-pulse method by inserting an extra 7r pulse into 
the light- pulse sequence. Instead of applying a 7r/2 pulse 
for the third pulse, the two atomic wavepackets are al- 
lowed to pass through each other and form a second loop 
[see Fig. |l0](b)]. Next a second tt pulse redirects the 
wavepackets to close the second loop, and a final 7r/2 
pulse interferes the wavepackets. This pulse sequence di- 
rectly produces a phase shift proportional to the the grav- 
ity gradient by essentially performing a coherent subtrac- 
tion of two spatially, and temporally, separated gravity 
measurements (one for each loop). Following the pulse 
rules presented in Section I, the double-loop phase shift 
is given by Acf> — <fi(ti) — 2<\>(t2) + 2(f)(t^) — ^(£4), where 
<p(U) is the phase of the i th Raman pulse at the position 
of the atomic wavepacket at time t^ . For a given fountain 
height the double loop interferometer is maximally sen- 
sitive when the interferometer spans one half of the total 
fountain time. Evaluating this phase, assuming a linear 
gravitational gradient g(z) = g + az along the sensitive 
axis, gives Acj> « 8kg aT 4 , where g is the gravitational 
component along the Raman wavevector at the atoms' 



initial position and a is the linear gradient. This formula 
is valid for an interferometer in the first half of a fountain 
of length 8T and for az' <C g D where z' is the total height 
of the fountain. 

One problem with this approach is that if the inter- 
ferometer spans the whole fountain time, i.e. is sym- 
metric about the fountain's peak, then the gravity signal 
from each loop will be identical, resulting in no phase 
shift. Thus, to obtain gradient sensitivity, the double- 
loop must be used only in the first half or the second half 
of the fountain, which severely constrains the possible 
instrument sensitivity. (This means the maximum inter- 
action time for our fountain is T ~ 39 ms, resulting in 
3000 times less sensitivity to gradients as compared with 
two single loop interferometer which are separated by 10 
m.) In addition to this sensitivity limit, it is difficult 
to make the gradiometer baseline arbitrarily large and 
there is no common-mode vibration rejection, since each 
of the two acceleration measurements which comprise the 
gravity gradient signal are made at differing times. 

A slight modification of the double-loop method al- 
lows full use of the fountain interaction time, resulting 
in a three-fold increase in sensitivity from a double-loop 
interferometer in an apparatus of equivalent size. This 
modified sequence uses an extra tt pulse ina7r/2 — 7r — 
7r — tt — ir/2 sequence, which creates a triple-loop in- 
terferometer as seen in Fig. O(c). The phase shift is 
Ac/) = - 2</)(h) + 2</>(t 3 ) - 2(j){U) + <t>{h), which 

reduces to: A<j> sa [17/3 + 4V2]kg aT 4 , again assuming 
a linear gradient and an interferometer that now spans 
the entire fountain time, (8 + 4y2)T. This interferometer 
may be symmetrically spaced about the fountain's apex, 
so that the maximal interrogation time is T = 46 ms for 
our fountain. Interference fringes from a triple loop in- 
terferometer are shown in Fig. [ll]. The 10 m single-loop, 
two chamber gradiometer still is 1000 times more sensi- 
tive due to its larger baseline and longer effective inter- 
rogation time. In addition, the triple-loop gradiometer, 
like the double-loop, does not provide vibration rejection, 
has a limited baseline, and is only slightly less sensitive 
to magnetic fields. 

2. Curvature Measurements 

Two multiple-loop gradiometers may be used in con- 
junction to measure the second moment of the gravity 
field, in a configuration similar to that used to measure 
the gradient with two single loop interferometers, as il- 
lustrated in Fig. H In this case, the triple-loop gradient 
phase shifts obtained from two spatially separated atomic 
ensembles are subtracted to obtain the second order cur- 
vature of the gravitational field. This device does have 
immunity to spurious vibrational noise because measure- 
ments are made simultaneously with respect to a common 
platform. Proof-of-principle data is shown in Fig. [ll]. 
Measurement of the second moment of the gravitational 
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field allows differentiation between massive, distant ob- 
jects and less massive, nearby objects that our gradiome- 
ter could not distinguish. 

A more direct and sensitive way to characterize the 
second order moment is to operate simultaneously three 
single-loop accelerometers spaced equidistantly along a 
single axis. The difference between gradients obtained 
by differencing the first and second accelerometer outputs 
from the second to third accelerometer outputs gives a 
second moment measurement. This device should have 
the same benefits as the previous device and would be 
significantly more sensitive. 

V. CONCLUSION 

Future sensitivity enhancements are likely. Improve- 
ments in the number of atoms trapped and optically pre- 
pared in the correct state would allow for higher, atom- 
shot noise limited SNRs. This might be accomplished 
using recently demonstrated atom trapping and cooling 
techniques pUJslfl . Multiple pulse (> 2Kk) techniques 
will provide straightforward access to large-area, high 
sensitivity configurations. 

It should be noted that the performance of a grav- 
ity gradiometer in a microgravity environment would be 
greatly enhanced due to the larger available interrogation 
time (A(/> ~ T 2 ) without need for an atomic fountain, and 
such a high sensitivity device would be ideal for tests of 
fundamental theories. 

Finally, the current sensitivity makes measurement of 
geophysical gradient signals possible. A portable abso- 
lute gradiometer would be useful for navigation, geodesy, 
terrain estimation, and oil and mineral exploration. 
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FIG. 1. A comparison of the transfer efficiency of a com- 
posite 7r pulse (solid traces) with a regular ir pulse (dotted 
traces) . The detuning is from the F = 4 m/ = 0->F = 3mj 
= transition. 

FIG. 2. Schematic illustration of the detection apparatus 
showing the Raman beam racetrack setup. 

FIG. 3. The decrease in the gradiometer SNR is shown due 
to a phase mismatch induced by the Earth's gravity gradient. 
The phase shift A(j> is the gradient phase shift from increasing 
the interferometer time T (A(f> txT 2 ). The solid line is a 
theory based on the noise from two mismatched sine waves in 
quadrature with photon shot noise from the detection system. 



FIG. 4. Gradiometer interference fringes with T = 157.5 
ms. Squares and circles represent the normalized upper and 
lower chamber fringes respectively, after compensating for 
contrast differences between the chambers. Triangles show 
the lower chamber signal residuals. 

FIG. 5. Schematic of the vibration isolation apparatus fea- 
turing elements for passive and active isolation. 

FIG. 6. Gradiometer interference fringes with T = 157.5 
ms. Circles represent the raw data points, and the solid lines 
is a sinusoidal least squares fit. (a) shows a fringe with the 
inertial reference platform, i.e. the cornercube, rigidly at- 
tached to the optical table, (b) shows an interference fringe 
taken with the cornercube mounted on the vibration control 
system. 

FIG. 7. Results of the reference platform shake test. No 
SNR reduction is seen when driving accelerations on the ref- 
erence platform at amplitudes of 10~ 2 g over the frequency 
band indicated. This amounts to 140 dB of vibration rejec- 
tion. The observed maximum SNR is slightly reduced from 
the data of Fig. || 

FIG. 8. Results of the platform tilt test. The floated op- 
tical table is tilted at a number of frequencies, and the gra- 
diometer is insensitive to a range of tilt amplitudes. The 
observed maximum SNR is slightly reduced from the data of 
Fig. | 

FIG. 9. Gravitational tidal signals as monitored by one ac- 
celerometer output over two days. Data are scattered points 
and the solid line is the tidal model with no free scaling pa- 
rameters. 

FIG. 10. Recoil diagrams for various interferometer pulse 
sequences. Dotted lines represent the F = 4 states and solid 
lines are the F = 3 state. The sensitivity to accelerations is 
proportional to the area enclosed in the recoil diagrams. The 
timing of n and 7r/2 pulses are shown with the vertical arrows, 
a) Single-loop accelerometer. b) Double-loop gradiometer. c) 
Triple-loop gradiometer. 

FIG. 11. Typical interference fringes from T = 44 ms triple 
loop interferometers in the upper and lower chamber. The 
solid lines are least squares fits. The phase difference be- 
tween the two fringes represents a measurement of the second 
moment of the gravitational field. 
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